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We consider the Anomalous Hall Effect (AHE) in thin disordered ferromagnetic films. Using a 
microscopic model of electrons in a random potential of identical impurities including spin-orbit 
coupling, we develop a general formulation for strong, finite range impurity scattering. Explicit 
calculations are done within a short range but strong impurity scattering to obtain AH conductivities 
for both the skew scattering and side jump mechanisms. We also evaluate quantum corrections due 
to interactions and weak localization effects. We show that for arbitrary strength of the impurity 
scattering, the electron-electron interaction correction to the AH conductivity vanishes exactly due 
to general symmetry reasons. On the other hand, we find that our explicit evaluation of the weak 
localization corrections within the strong, short range impurity scattering model can explain the 
experimentally observed logarithmic temperature dependences in disordered ferromagnetic Fe films. 
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I. INTRODUCTION 

It has been recognized since the 1950's [H that a Hall 
effect can exist in ferromagnetic metals even in the ab- 
sence of an external magnetic field, hence the name 
Anomalous Hall Effect (AHE). There are several differ- 
ent mechanisms that might be responsible for the AHE 
observed in thin ferromagnetic films, namely the skew 
scattering Q and side jump mechanisms [1] as well as 
Berry phase contributions All such mechanisms de- 
pend on spin-orbit interaction induced by the impurities 
and on the spontaneous magnetization in a ferromag- 
net which breaks the time reversal invariance and there- 
fore gives rise to the AHE. For a disordered ferromag- 
netic film, AH conductivity due to the skew scattering 
and side jump mechanisms have been theoretically con- 
sidered using a variety of methods within weak, short 
range impurity scattering 0, H, 0, BUl- However, a sys- 
tematic calculation, starting from a microscopic Hamilto- 
nian, of the longitudinal as well as the AH conductivities 
for different mechanisms for strong impurity scattering 
has been lacking. Recently, the effect of strong, short 
range impurity scattering on the longitudinal and Hall 
conductivities were considered for skew scattering as well 
as side jump mechanisms [10], but quantum corrections, 
namely electron-electron (e-e) interaction corrections [1 1] 
or weak localization (WL) effects [l^l , were not included. 

Earlier experiments [13| have shown logarithmic tem- 
perature dependences of the longitudinal as well as Hall 
resistances highlighting the importance of such quantum 
corrections. However, the results were consistent with. 
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and were interpreted as, vanishing interaction contri- 
butions to the AH conductivity, obtained theoretically 
within a weak impurity scattering model fl| and the ab- 
sence of any weak localization effects. Recent experi- 
ments on the other hand clearly show non- vanishing con- 
tribution to the total quantum correction to the AH con- 
ductivity , which can arise in principle either from an 
interaction correction due to strong impurity scattering, 
or from a weak localization effect, or from a combination 
of both. It has been commonly believed that weak lo- 
calization effects in ferromagnetic films would be cut off 
by the presence of large internal magnetic field among 
others, which suggests that the interaction corrections to 
the AH conductivity need to be revisited for strong im- 
purity scattering as a source of difference between the 
two experiments. 

In this paper we systematically develop a general for- 
mulation for the AHE for strong, finite range impurity 
scatterings starting from a microscopic model of electrons 
in a random potential of impurities including spin-orbit 
coupling. This generalizes an earlier work Q which con- 
sidered weak, short range impurity scattering only and 
did not include quantum corrections. We show on very 
general symmetry grounds that quantum correction to 
the AH conductivity due to (e-e) interaction effects van- 
ish exactly, which shows that the previous weak scatter- 
ing results ^ remain valid for arbitrary strengths of the 
impurity scattering. This forces us to consider the weak 
localization effects Q as the only remaining source of the 
logarithmic temperature dependence in the above exper- 
iments despite the presence of large internal magnetic 
fields and spin-orbit scatterings in these ferromagnetic 
films. As we show below, the temperature independent 
cutoff of the weak localization effects in strongly disor- 
dered systems can be ineffective at higher temperatures 
if a temperature dependent contribution dominates the 
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phase relaxation rate. It turns out that while the contri- 
bution from the e-e interaction to the phase relaxation 
rate is indeed too small for WL effects to be observed, 
a much larger contribution is obtained from scattering 
off spin waves jl^ , which should allow the observation of 
the WL effects within a reasonable temperature range. 
We find that the effects of strong impurity scatterings on 
the WL effects can be evaluated to obtain a very sim- 
ple result, namely that the ratio of the WL corrections 
to the AH to the longitudinal conductivity can be writ- 
ten simply in terms of the eigenvalues of the impurity 
averaged particle-hole scattering amplitude for zero mo- 
mentum transfer. This result, taken together with con- 
tributions to the AH conductivity from both the skew 
scattering and side jump mechanisms calculated within 
the same microscopic model, can explain both the ear- 
lier as well as the recent experiments on the disorder and 
temperature dependences of the AH conductivities of ul- 
trathin Fe films [ij] mentioned above. This last result 
has been reported without details in combination with 
the recent experiment in a short letter p^ . 

The paper is organized in the following way: A micro- 
scopic model Hamiltonian is introduced in Section H, and 
a general formulation in two dimensions for strong, finite 
range impurity scatterings is developed in Section HI. 
Section IV reviews the results on the conductivity tensor 
in the absence of interactions. In Sections V and VI we 
consider the e-e interaction corrections and the weak lo- 
calization corrections, respectively, to both longitudinal 
and AH conductivities within the general strong, finite 
range impurity scattering formulation. We then consider 
the special case of a short range, but still strong, impu- 
rity scattering model in Section VII. In Section VIII we 
collect all the results and compare them with recent ex- 
periments. Section IX summarizes the paper. For the 
sake of completeness, we include models of small and 
large angle scatterings in the Appendix. 



II. HAMILTONIAN 

The single particle Hamiltonian of a conduction elec- 
tron in a ferromagnetic disordered metal, including 
spin-orbit interaction induced by the disorder potential 
Vdis(r), is given in its simplest form by (throughout the 
paper we use units with h = = 1) 



Hi - 



(47r) 



{VVd^s X V)], 



(2.1) 



where Ac = ^ is the Compton wavelength of the elec- 
tron, and M is the Zeeman energy splitting caused by 
the ferromagnetic polarization. Here 7?i is a 2 x 2 ma- 
trix in spin space with ct, a' =| , | being spin indices and 
T is the vector of Pauli matrices. The above model is 
only a crude approximation of the bandstructure of Fe, 



which has been determined by several authors (see e.g. 
ref [13 )• We model the energy band crossing the Fermi 
surface by a single isotropic band. As will be discussed 
below, the quantum corrections to the conductivity ex- 
hibit certain qualitative features, which do not depend 
sensitively on the details of the band structure. The dis- 
ordered potential in (2.1) will be modelled as randomly 
placed identical impurities, Vdis(r) = V{v — Rj). We 
will later average over the impurity positions Rj . 

The matrix elements of Hi in the plane wave (or Bloch 
state) representation are given by 



(k/cr/|Fi|kcr) = j d 



'^re-'^'-^Hie-'^" 



' 2m 
V;o(k/cr/;kcr) 



(2.2) 



where F(k — k/) is the Fourier transform of the single 
impurity potential, and the spin-orbit interaction part is 
given by 



V;o(k/cr/;kcr) 



A? 



(47r)2 



^F(k-k/)eW''-''')-^^l 



X T^^i ■ (k X k/) 
Here we have used 



(2.3) 



y d^r exp(-ik/ • r) i^Vdis x V) exp(-«k • r) 



d^r 



^i(k-k/-q)-r 



(27r)2 
-iV{k-k/){k X k/) 



(-*q)V(q) X (zk) 



(2.4) 



The many-body Hamiltonian is given in terms of electron 
creation and annihilation operators c^^^jCucr as 

H = ^ (£k - A/a)c+^Ck, 

kcr 

+ J2 5]^(k-k/)e*(''-'''>^^- 

kcT,k/cr/ j 

X {Sera/ - igsoTcrcr/ ■ {kxk')}c^,^^Ckcr (2.5) 

where we have defined a dimensionless spin-orbit cou- 
pling constant = (^^rp-, k = k/|k|. Note: An es- 
timate of the spin-orbit coupling constant gso, using a 
typical Fermi wave number kp , shows that it is rather 
small, of order 10^"*. However, in transition metal com- 
pounds the couplingis substantially enhanced by inter- 
band mixing effects [SI] , so that the renormalized coupling 
constant gso is of order unity: gso ~ CsoEso/ ^Ed, where 
Eso ~ O.leV is a measure for the atomic spin-orbit en- 
ergy, AEd ~ 0.5eV is a typical energy splitting of d- 
bands, and the constant Cso ~ 5. In the following we 
will replace gso by the phenomenological spin-dependent 
parameter g^-. 
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III. IMPURITY SCATTERING: GENERAL 
FORMULATION 

In this section, we will develop a general formulation 
for strong, finite range impurity scattering in two dimen- 
sions using standard field theory techniques at finite tem- 
perature jlTj. For simplicity, we will need to make ap- 
proximations for short range impurity scattering later. 
However, keeping the formulation general as long as pos- 
sible will allow us e.g. to check if the anisotropic scatter- 
ing can have a large impact on our final results. 

The repeated scattering of an electron off a single im- 
purity may be described symbolically in terms of the scat- 
tering amphtude fka,k/ai as 

f + VGV + VGVGV + .... (3.1) 

where G is the single particle Green's function 

Gka{it^n) = [iUJn - Ska - Sfco-(w„)]"\ (3.2) 

with the single particle self energy T,ka-{ii^n)- Here w„ = 
7rr(2n + 1) is the fermion Matsubara frequency with T 
being the temperature and A^cr is the density of states 
at the Fermi level of spin species a. (We use units of 
temperature such that Boltzmann's constant is equal to 
unity). V is the bare interaction with one impurity at 
R = and includes the spin-orbit scattering 

VkMa = V{k - k')[l - ig,Tl„(k X fc')], (3.3) 

where we have used the fact that V is diagonal in spin 
space. In the case of finite range, or even long-range cor- 
related scattering potentials, we may still use the model 
of individual impurities or scattering centers, but now of 
finite spatial extension. This is reasonable as long as the 
scattering centers do not overlap too much. If they over- 
lap, a more statistical description in terms of correlators 
of the impurity potential should be used. Within our 
model, the nonlocal character of scattering is described 
in terms of the momentum dependence of the Fourier 
Transform of the potential of a single impurity (assum- 
ing only one type of impurity) l^(k — k/), which for an 
isotropic system depends only on the angle d between k 
and k', V = V{&) = Vy-O). In 2d we may expand V 
in terms of eigenfunctions Xm(A) = e™"^, where is the 
polar angle of vector k, fc = k/|k|. Adding the skew 
scattering potential we may write 

Vk,k,a = Vrn.Xr,r(k)xUk') (3.4) 

rn 

where Vmcr is a sum of the normal and skew scattering 
parts 

Vrr..^V„r+VZ- (3.5) 

Time reversal invariance and rotation symmetry in the 
case of potential scattering implies 

= (O* = C'- (3-6) 



Equation (3.3) then yields 

VZ = Ig.r^^ (Cli - Cii ) (3.7) 

A. Scattering amplitude 

For V diagonal in spin space, the scattering amplitude 
fka,kici = Scr,a' fk,k/o- obcys the integral equation 

fLk/a = VkMc + y^,Gkicr{iuJn)Vk,kicrfk,,kla- 
ki 

= Vk^T '^isTTN^{Vk,kiafkiM<y^ki, (3.8) 

where s = sign{ujn) and (. . denotes averaging over 
the direction of wavevector ki. Defining dimensionless 
potential Vma = T^N^Vma and the dimensionless scat- 
tering amplitude fka,kia/ = T^Na-fka,kia/ and expanding 

fkaMla = Em fmaXnAk)x'^{k') , WC find 

= ■ (3.9) 

For notational simplicity, we will always use a bar on 
a symbol to represent the corresponding dimensionless 
quantity. 

B. Single particle relaxation rate 

The single particle relaxation rate is given by the 
imaginary part of the self energy 

= -sn,mplm{f^^,,^) = -^la (3.10) 

where is a dimensionless parameter characteriz- 
ing the scattering strength, "/a- = -sY.m^'^Uma) = 

J2m i_|_'y2 J and Ncr is the density of states at the Fermi 

energy of spin species a. Note that Vma are all real. 

C. Particle- hole propagator 

The particle-hole propagator T^k/iQ] i^n ~ i^m) is 
an important ingredient of vertex corrections of any kind. 
Here k + q/2, k — q/2 are the initial , kf + q/2, kl — q/2 
the final momenta and e„,e„ — are the Matsub- 
ara frequencies of the particle and the hole line, respec- 
tively. In terms of the particle-hole scattering amplitude 
tk,kr{<l','i'£m'i^m), F Satisfies the following Bethe-Salpeter 
equation (we have defined dimensionless quantities F, t 
by multiplying both with a factor {27rN„T„)) 
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+ {2'KNaT„)~'^ ^ tkki {q; ie-n, i^m)Gk^+q/2,a{i^n) 
fei 

X Gkj^_q/2,a{ien - i^m)^kik/{q;ien,i^m) (3-11) 

The (dimensionless) impurity averaged particle-hole scat- 
tering amplitude t (we consider only the case of equal spin 
of particle and hole) is given in terms of the (dimension- 
less) scattering amplitudes / by the equation 

Tssl j O \ '^'^c'^^imp 7s ( ' \ 

t^kk/yq^^^nj^^l'm) — ./fc+g/2,<7;fc/+q/2CTV*£raJ 

^ fkl-q/2,a;k-q/2,cri'^^n — i^m)- (3-12) 

We will later need the limit of small q, q <C kp, of this 

expression, 

tl%iq; *e„, = t^^ijl = 0) + A PM- (3.13) 

It is useful to represent the operator ikk/{q = 0) in terms 
of its eigenvalues Am- Assuming isotropic band struc- 
ture, the eigenfunctions Xm{k) = ex.p{invp) are those of 
the angular momentum operator component L^- The 
eigenvalue equation is 

{tkk,iq = 0)Xm(fc'))fe' = ^mXmik). (3.14) 

The operator i^^, {q = 0) may be represented as 



t 



*fcfc!(9 = o) = [i+-k{q = o)r (3.15) 
Ig 



"kki Vy — — I "fc/fc 

In general, using the definitions 

^imp Js 



(7riV^)2 



fk,k'<Tfk',ka — {^'l^^crTa) tff^k'ir 

(3.16) 



we have 



^m<7 ~ 7(7 fmlafml-m,a- (3-17) 



We will consider A ikk/{q) for the special case of strong 
short range impurity scatterings in section VII-A. 

The energy integral over the product of Green's func- 
tions in the integral equation for Tkki may be done first, 
after expanding the G's in Q^n and q, 



J 



dei Gkj^+q/2,a{'i^n)Gki-q/2,cr{i^n — i^m) 

= 27rr [1 + iTiiilrn - q • J 

- r\<i-Vkj'] (3.18) 



with e„ > and e„ — flm < 0, where q- Vj, = qvpiq- 

k). Expanding Tkki and tkk/ in terms of eigenfunctions 
Xm{k), Tkk/ = J2m^rnmiXm{k)x'^,{k') and using i+- = 
Xm one obtains (s' = —s) 

^m'ml = ^m^mm' + Ato{[1 - T{\fln\ + Doq^)]Tmrn) 



- ^VFqTs[T'^\^,X*l{q)+^"r^'+l,m,Xl{q)] 



- ^{vFqTr[Tt^'_2,m,X2{q) 

+ ^Z'+2,mMq)]} 

For m = m' ^ the solution is 
A. 



(3.19) 



r — 

■1- mm — 



l-A, 



+ 0{q) = Xm + 0{q), (3.20) 



where we have defined Xm = Am/(1 — Am). The A^ 
(and therefore A^) are complex valued and depend on 
the spin projection a. Using conventional notation, we 
will denote the real and imaginary parts of A^ by A^ and 
A^, respectively, and similarly the real and imaginary 
parts of Xm by A^ and A^, respectively. 

The case m = needs special consideration, because 
particle number conservation causes Too to have a pole 
in the limit f2„, g ^ 0, here expressed by Aq = 1. Solving 
the above equation for f oo in lowest order in q, one finds 



Too 



1/7 



(3.21) 



where the renormalized diffusion constant is defined as 

D = Do{l + A'l), Do = ^v%T 

A'l = ReAi = i(Ai-t-A_i); (3.22) 

This is found by solving the following equations for small 
VFqT{s' = —s) 

f-' = l + [l-T{\nm\+Daq^)no 

- -jVFqTs[T'':^^^-,x*l{q)+T\''^^Xl{q)] 

f%, = x.,{f^^i,-'-VFqTsfi^;,xim 
n% = xr{f{%-'-VFqTsfi^;ox-im (3.23) 

Substituting f j?{ q into the equation for f qq one finds: 

n^{\nm\ + Doq^[l + i(Ai + A-i)]} = - (3.24) 

Z T 

The leading singular dependence on ki is obtained from: 

fs:±i = [i-r(|0mi+i)o9')]fg:ii 

- \vFqTs[f^X^,xl{q)+ff!±iXim 

f!:ti = -l^X^iVpqTsT'^'I^Xliq) 

r?:( = Ai-^Aiz;fgTsf-;xt(?) (3.25) 

The complete particle-hole propagator in the regime 
vpqr < 1 \s given by 



fs 1 Iklk, 

i-kki = --, 



T \nm\+Dq^ 



E ^mXm{k)Xmik') (3.26) 
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with 



7fe = 1 - -^vpqTS ^ XmXm(k)X'>'n*{q) 

m=±l 

= 1 - ^VpTS ^ XmXm{k)q-m (3.27) 



m=±l 



and 



7fc = 1 - l^'VpTS XmX*m(k)qm- (3.28) 



m— ±1 



The vertex corrections of the density and current ver- 
tices jka and jka (for the incoming and outgoing current) 
are obtained by 



Tk{q) = 1 + (f = 1 



l/r 



7fe (3.29) 



and 



3ka{q) = Vka + {Vkla^klk)kl 



m=±l 

l/r ^ 



+ {Vklalkl) 



kr. 



\^m\+Dq^ 
3ka{q) = Vka + {Vkla^kkl)kl 



Ik 



Vka + {'k){Vk,aX*m{k'))kl 



m=±l 



+ {Vktalk)kr 



l/r 



7fe- 



(3.30) 



Note that jka 7^ ( jka)* 1 as the eigenvalues are in 
general complex valued. Using 

X-i{k)xl{k') + Xi(k)x-i{k') = 2{k-k') 
X-i{k)xUk') - Xi{k)X-i{k') = 2i(kxk') (3.31) 



and 



k'aik-k')) = ^kai {K(kxk')) = ~{ea X k) (3.32) 



and defining jka = jka{q = 0), jka = jka{q = 0), we have 

ha = wf[(1+ A'Jfc„ + >/l{eaXk)^] 

jka = VF[il+ K)ka- X'lieaXk),]. (3.33) 

More explicitly, for a = x,y, the incoming and outgoing 
current vertices j and j have the forms 

jkx = w_f[(1+ X[)kx+ X'(ky] 

= Ivf[{1+ Xl)k+ + {1+ Ai)fc_] 

jky = vf[{^+ K)ky- A'/fcj.] 

= -^\vf[{1+ K)k+-{1+ Ai)fc_] 



jkx = vf[{1 + X'-{)kx - X'Iky] 

= Ivf[{1+ Xi)k+ + {l + Xl)k-] 

jky = Vf[{1+ X[)ky+ X'[k:r] 

= -i^VF[{l + Xi)k+ - (1 + Xl)k-] (3.34) 
where we have defined k± = k^i iky. 

D. Particle-particle propagator 

The integral equation for the particle-particle propaga- 
tor or Cooporon reads (again multiplying the Cooperon 
C and the particle-particle scattering amplitude by the 
factor 2'KN„Ta to define dimensionless Cooperon C and 
dimensionless particle-particle scattering amplitude tP) 

+ (27rA^crrcr)~^ y^t|fc^ {Q;ien,i^rn)Gki,cr{ien) 

kl 

X GQ-ki,(T{ien-i^m.)Ckiki{q;ien,i^m) (3.35) 



,p,ss 
''k,kla 



I'imp 7s 7s' 

I J k,k' a J —k, — k\a 



{2nN,T,)-'-f-'rk,k'J-k,-k',a (3-36) 



t^^a = '^^N^r.tlii = E e/'x™(fc)X™(^') (3.37) 



fP,Ss' ^ -1 ns ns' 

^mcr la / ^ J mfa J m—ml ,<7 



(3.38) 



If rotation invariancc or time reversal invariance is bro- 
ken, F^f = 7-17^ ^ 1, where 7? = f^,af^m,,a- 

The energy integral over the product of Green's func- 
tions in the integral equation for Ckk/ may be done first, 
after expanding the G's in flrn and Q , 



/ 



dSl Gki,a{i^n)GQ-ki,a{'i£n - i^m) 

= 2ttt[1 + iT{iQ,m - Q • J 

- r\Q-v,J% (3.39) 



with e„ > 0, e„ — flm < 0, where Q • Vj, = Qvf{Q ■ 
k). Expanding Ckki and in terms of cigcnfimc- 

tions Xm{k), Ckki = J2m^mmrXm{k)Xm'{k') and denot- 
ing = A^ one obtains (s' = —s) 

Cmml = X^{5„irn' + [1 ^ r(|0„| + DaQ^)\Craml 
X — ^ _ ^ 

- :^WFQr[C™_i,„,x*((3) + Cm+i 

,m/Xl m 

- ^{vFQTf[Cm-2,m.lX2{Q) 



Cm+2,miX2 



(3.40) 
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The m = to' = component of Cmmi obeys the equation where Lap = X] -^^p ^ sum of the different relevant 

n 

diagrams d„. We wiU take the current to be along x di- 
rection, so the longitudinal conductivity will correspond 
to a — P = X while the (anomalous) Hall conductivity 
(3 41) ^^^^ given by the off-diagonal part a = x, (3 = y. Note 

that Cyx — ^xy- 



[(T^°)-i + |r!„| + i?og']Coo = r-i 

o(O') 



where {t^°) ^ is the phase relaxation rate contributed by 
spin-orbit interaction processes: 



-M(Ag)-i-i]. 



(3.42) 



Using 



C±i,o = ^li{C±ifl - QrCcoXii®} (3.43) 



the Cooperon is found as 

Gkkl — 



p~rs 

7fc7fe, 



1 — Am 



(3.44) 



with 



A. Skew scattering contribution 

The skew scattering contribution to the conductivity 
tensor aap in lowest order in l/epT is given by the 
bubble diagram dressed by vertex corrections given by 
the correlation function 

Lap = T^^Gka{ien)Gka{ien - i^lrn)vkajkp (4-2) 

The energy integration over GG is nonzero only if the 
poles are on opposite sides of the real axis, requiring < 
e„ < flm (we assume flm > ) and yields 2t: N^t^, and 
the summation on e„ gives r2„i/(27rr). Substituting j^^ 
from Eq. (3.34) into the Kubo formula, the conductivity 
tensor follows as 



p 
Tk 



1-\vfQt E KiXm.(k)X*yn{Q) 



m— it 1 



ss _ 1 2 1 + A'l A'/ 



(4.3) 



- E KnXAk){Q,-^kiX'm[k')) (3.45) Defining the tensor of diffusion coefficients as 



m=±l 



and 



Tk 



] = A^;,x™(fc)(q-v„x™(fc')>-(3.46) 



m=±l 



Here the diffusion coefficient is in general different 
from the one in the p-h channel, 

i^P = 7^o[l + ^(A?-f A^)]^ A (3.47) 

the difference being proportional to the spin-orbit cou- 
pling Qa- 



D 



2' ^ 

Dly = Dl^[X'i/{l + \',)]^~Dl^ (4.4) 



where 



r.(l + A'i) 



(4.5) 



is the momentum relaxation time, we may write 

<%^Y.^^DU- (4-6) 



IV. CONDUCTIVITY TENSOR IN THE 
ABSENCE OF INTERACTION 

As mentioned before, there are three mechanisms con- 
tributing to the anomalous Hall conductivity, namely the 
skew scattering, the side jump and the Berry phase mech- 
anisms. In this section we will write down the generic for- 
mulations for evaluating these contributions within dia- 
grammatic perturbation theory. The contribution to the 
conductivity Gap will be given in terms of a correlation 
function Lap, defined as [l3| 



Ca/s = lim— 



Sl-»0 



(4.1) 



From the definition Am = Am/(1 — Am) we obtain the 
following identities: 

1 + Ai = 1 + a;- ^"^1 



11 -All 



1-Ai 



B. Side-jump contribution 

The side-jump contribution has been first calculated 
by Berger [3|]. It arises because the trajectory of a wave 
packet scattered by an impurity is shifted sidewise due 
to the spin-orbit interaction ("side-jump"). This effect 
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FIG. 1: Diagrams for side jump contributions. Solid lines 
are impurity averaged Green's functions. Shaded triangles 
with dashed lines represent impurity scattering amplitudes 
while the dotted line from a vertex denotes spin-orbit term 
in the velocity operator. The shaded vertex represents vertex 
corrections to the current density operator. 



[101 . Analogous contributions for thin film ferromagnets 
have not been obtained yet. Such contributions depend 
on the details of the band structure and is beyond the 
scope of the present work. On the other hand, the fo- 
cus of the current work is on the disorder and temper- 
ature dependence of the AH conductivity in which the 
Berry contributions are qualitatively similar to the side 
jump contributions (both arise from an additional veloc- 
ity term due to spin-orbit interactions). Therefore, the 
effects of Berry contributions can be included in a phe- 
nomenological way, while comparing with experiments, 
by considering a larger side jump contribution to the to- 
tal AH conductivity. 



may be calculated in a straightforward way [l8| by ob- 
serving that the side-jump leads to an additional term 
in the particle velocity due to the spin-orbit interaction. 
Indeed, the quantum mechanical velocity obtained from 
the Heisenberg equation of motion for the position oper- 
ator has two terms. 



^ T 17 1 P 

V =— r = —i[r,Hi\ = — 



1 



;iTXVVd^s)■■ (4.. 



dt '''''' m Am?c^ 

The Bloch states matrix elements of v are given by 



(k/cr/|v|kCT) = —5^\^,5„o 

m 



9o 



2ni£p 



E^(k-k/) 

(4.9) 



For strong impurity scattering, there are six diagrams 
that contribute to the current correlation function, four 
of type (a) and two of type (b), shown in Figure 1. For 
example, contributions from diagrams of Figure 1 (a) and 
(b) give 



k-k' 



kk' 



k-k' 



kk' 



^ fk'kifkikJkiy 



2m 
(4.10) 



These were evaluated within the short range strong im- 
purity scattering model in Ref. [l3|. We will later use 
the results reported there. 



C. Berry phase contribution 

In general. Berry phase contributions can arise when 
there is an anomalous velocity term, as in the case of 
the side jump contribution given by Eq (4.8). In princi- 
ple, such terms can also arise in the presence of a peri- 
odic potential and spin-orbit interaction leading to finite 
Berry curvatures It has been found that the intrin- 
sic Berry curvature contributions to the AH conductivity 
for bulk ferromagnetic metals can be large in magnitude 



V. INTERACTION CORRECTIONS TO THE 
CONDUCTIVITY 

The e-e interaction corrections to conductivity will be 
calculated in first order in the screened Coulomb interac- 
tion. It may therefore be represented as an integral over 
a kernel K(q,iuji) multiplied by the screened Coulomb 
interaction Vc{q, ii-^i). 



6a' = TJ2J dq^Kiq, iiUi)V,iq,iu;i). 



(5.1) 



Gauge invariance requires that Sa should be invari- 
ant against an energy shift of the interaction potential, 
V{r) ^V{r)+C, which only leads to a constant term in 
the total Hamiltonian. In Fourier space, the transforma- 
tion is V{q) — *-y(q) -I- C(5(q), which requires the kernel 
to vanish in the limit q ^0 (Even more general, 

since V{q) is an electric potential, a gauge transformation 
of the above form, but with arbitrary time dependence 
C = C{t) does not change the physical fields.) We will 
see below that this gauge invariance, together with an ad- 
ditional mirror symmetry, will impose a strong constraint 
on the interaction corrections to the Hall conductivity. 



A. Coulomb interaction renormalized by diffusion 

The Coulomb interaction Vc((7, cj/) is renormalized 
by diffusion processes. The bare screened interaction is 
given by 

Vc{q,iu;i) = VB{q)/[l + VB{qMq,iiUi)], (5.2) 

where Vsiq) — 47re^/g^ in 3d and Vsiq) = 27re^/g in 
2d, and the polarization function is given by [l^ 



H(g, iuji) 



dn Dq^ 
dfl \iOi\ + Dq^ ' 

In 2d one therefore finds 

27re2 \u;i\+Dq^ 



Vc{q,iu!i) 



q \uji\ + Dq^ + DqK2 ^ d^i 



(5.3) 
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(a) 



(b) 



FIG. 2: Diagrams for interaction corrections. Solid lines 
are impurity averaged green's functions, wavy lines denote 
screened coulomb interactions and dashed lines denote diffu- 
sion poles. There are two diagrams of type (a) and two of 
type (b). 



in the limit w; = 0, g — > 0. Note that in a ferromagnet, an 
additional effective electron-electron interaction arises by 
exchange of spin- wave excitations. We do not consider 
this interaction here because it is small, of order (J/ef)^ 
where J is the exchange energy (see section VI-B). 



B. Singular contributions for skew scattering 

The diagrams for the correlation functions Lap de- 
fined in (4.1) can have up to three diffusion poles [2]| . 
The gauge invariance argument presented above suggests 
that the relevant contributions to K{q,iuji) should have 
a factor of q^, which cancels one of the diffusion poles. 
Therefore only diagrams with three diffusion poles shown 
in Figure 2 contribute. For example, contribution from 
diagram (a) of Figure 2 is given by 

Ll^P = -^E^E E G2(e„)Gfc_,(e-^) 

Cn 1^1 k,k',q 

X Gk'~q{e - w)G'fc'(e„)Gfe'(e„ - n)V{q,uJi) 
X [e(e)e(e-r!)e(c^-e)T+-(q,c^) 
X T,-+(-q,-c.)r+^(q,c.-r!) 



X T+-(-q,-c.)r+J(-q,^ + r!)] 

X VkaVk'13- 

Using only the singular parts 



(5.5) 



r+-(r, o) - 7fc(q)7fc'(q) 
rfefet(q,w) = r+j;(-q,-w); 



(5.6) 



and 



r,7+(q,^) 



and defining 



7fc(q) 

\uj\+Dq^ 

7fc(-qJ 

H+Dq-^ 



V{q,iOi) 



i\u;i\+Dq^)^{\ui-nr, 



Dq^ 



(5.7) 



(5.8) 



one gets 

Ll'p - Y.i-2mNoT^rY.^T ^ (c., - a„) 

X {vkalk (q)7fe (q)Cfe (q)) {vk'f3ik' (q)7fc' (q)C/c' (q)) k' 
+ k;|(t'/ca7fc(-q)7fc(-q)4*(q)>fe 

X {vk'pjk' i-qhk' (-q)Cfc' (-q))fe'] 

X -^Pq(c^,,rj„) (5.9) 

where we have expanded the Green's functions for small 
q and defined the factor 



= 1 - 2ir(q • Vk). 



(5.10) 



Note that 7fc(— q, —^) = 7fc(q, ^) The leading terms in q 
are the linear in q terms in the products 77^: 



7/c(±q)7fc(±q)Cfe(q) = l T 2iT(q • Vk) 



(5.11) 



m=±l 



The A's combine to X'^ — A'_„, which may be pulled in 
front of the m-summation. Observe that 

VF E Xm(^)g-m = 2(q • Vk). (5.12) 
m=±l 

Therefore quite generally, 

(wfcx7fe(q)7fc(q)C/c(q))fc = -ivprq^il + X[). (5.13) 



C. Corrections to longitudinal conductivity within 
sltew scattering model 

For contributions from diagram (a) of Figure 2 to the 
longitudinal conductivity, each of the two angular aver- 
ages (in each term) in Eq. (5.9) with a = (3 = x gives 
a factor proportional to qx (see Eq. (5.13)), the product 
yielding q^. Diagram (b) also has the same combination. 
This yields, for the four diagrams (a), (a'), (b) and (b') 
the total contribution (i^^ ^ ]^2a\^ i2b^ ^ ]^2b^y 

r 2a+2a' +2b+2b' ^ S^/o^AT _2^2/ 2\2/^ , 7l \2 

Lxx = ^ Z^(27rAr„r ) [vpT) (l + Aj 



(5.14) 



where we have defined 



■^{q,^m) = TY,^lP{-^U^m)-'D{-OJl-^m,^m)] 
W!>0 



(5.15) 
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The sum over q converted to an integral yields 

E'^^*('^'^™)-il^^(l + l-^)' (5-16) 
q 

where k = 27re^ N^r is the screening length. 

The exchange interaction correction to the longitudinal 
conductivity is then given by 

S<^T. = ^L,. = -^^^'^, (5.17) 

where we used D„ = D()a{\ + X'i„). Note that the cor- 
rection 5<Txx is independent of scattering strength. 



D. Corrections to Hall conductivity within skew 
scattering model 

For a = X, f3 = y, the two angular averages in (5.9) 
are proportional to Qx and qy, respectively, so that the 
angular g-integral yields zero. This is true for all four 
diagrams (a), (a'), (b) and b'. Thus the total correction 
to the Hall conductivity L^y within the skew scattering 
model is zero. Note that the results are true for arbi- 
trary strength as well as finite range and anisotropy of 
the impurity scattering. 

Note that the result that the angular average (5.13) 
is proportional to qx is a special consequence of the fact 
that (5.9) contains the combination ^klk- This partic- 
ular combination is proportional to q • Vk as shown in 
(5.12), which results in (5.13). This is true for the class 
of diagrams considered here. This leads to the obvious 
question if there are other diagrams where the angular 
average is over a different combination of 7^ 's leading to 
a non-zero contribution to Lxy It turns out that indeed 
there are such terms with less than three diffusion poles, 
but that there is a deeper reason why the total inter- 
action correction to the Hall conductivity must always 
vanish in the first order in Coulomb interaction. In this 
case, the interaction correction has the form (5.1) and the 
kernel must be proportional to q^ as mentioned before. 
In addition, we have the following symmetry properties 
for the Hall conductivity with respect to a sign change of 
the magnetization (magnetic field) and a mirror reflec- 
tion from the yz-plane x —x (or from the sz-plane 
y —y) which follow from the invariance of the Hamil- 
tonian under a simultaneous transformation B ^ —B 
and X — > —X (or y — > —y): 

<Jxy{B) = -axyi-B) 
axy{B]x) ^ axy{-B\-x) ^ -(Jxy{B;-x) (5.18) 

which means that the Kernel must be proportional to 
qxqy to preserve the mirror symmetry. Thus, even though 
individual diagrams do contribute, the total sum of all 
diagrams of a given class must cancel to yield vanish- 
ing contribution to the Hall conductivity. Note that the 



above argument remains valid for the side jump contri- 
butions as well. Therefore we have, quite generally, 

5aiy = 0. (5.19) 

This generalizes the results of LW where this result was 
first obtained within a skew scattering model with short 
range and weak impurity scattering. 

Note that the above arguments do not imply that 
the weak localization correction to the Hall conductiv- 
ity must also vanish, because the WL contributions do 
not have the form Eq, (5.1) and the gauge invariance 
arguments do not apply. 

E. Corrections to conductivity within side jump 
model 

We have already argued that the e-e interaction correc- 
tions to the Hall conductivity due to side jump scatter- 
ing must vanish on very general symmetry grounds. The 
corresponding corrections to the longitudinal conductiv- 
ity are of course finite. However, these contributions are 
proportional to the spin-orbit coupling, and therefore are 
much smaller than the corrections due to normal scat- 
tering obtained above. We will therefore neglect such 
contributions. 



F. Hartree terms 

Eq. (5.17) should be corrected by including diagrams 
of the Hartree type. This leads to the total interaction 
correction in 2d [ll| 

^'^- = -^(l"i^'^)l^y' (5-20) 

where 

= 8(H-F/2)ln(l + F/2)/F-4 (5.21) 

and 

F= , \, / ^v{q^2kF&\ne/2). (5.22) 
VKQ = 0) J 27r 

As we will discuss later, experiments suggest an approx- 
imate cancellation between the exchange and Hartree 
terms, which will imply that the quantity 

hxx = (1 - ^F) (5.23) 

can be very small. 

VI. WEAK LOCALIZATION CORRECTION TO 
CONDUCTIVITY 

As pointed out before, the weak localization contri- 
butions can not be written as an integral over a kernel. 
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FIG. 3: Diagrams for weak localization corrections. Solid 
lines are impurity averaged green's functions, broken lines are 
impurity scattering amplitudes. Shaded cross is the cooperon 
and shaded vertices axe vertex corrections to the current den- 
sity operator. There are two diagrams of type (b) and four 
diagrams of type (c) 



as in (5.1) for the Coulomb interaction. Therefore, al- 
though the mirror symmetry is still preserved, the total 
contribution to the Hall conductivity need not be zero. 



A. Cooperon contributions 

The weak localization correction to the current- current 
correlator is obtained from diagrams shown in Figure 3, 
with one Cooperon propagator connecting the upper and 
lower line of the conductivity bubble. The frequency ar- 
guments of the upper (particle) line and the lower (hole) 
line have opposite signs. The current vertices are dressed. 
For example, the contribution of diagram (a) of Figure 3 
to the current correlation function is 



r 3a 



<y e-n k,k/,Q 

X {2T:N„Ta)~^Ckk,{Q\ien,i^m) (6.1) 



Here the momentum Q = k -I- k' can be taken to be 
small, as for Q ^ the Cooperon is strongly peaked. 
Consequently one may take k' k, — k in the arguments 
of the Green's functions and of the current vertex, i.e. 
^ -jki3- Then 



T 3a 



(r!™/27r) ^(47rA^,T^)(27riV,T.)-i 



X {jLrk0)kJ2^kMQ) 

Q 

The Cooperon contribution is given by 



(6.2) 



Q 

'"^^ QdQ 



1/r 



JO 2^ |n„,|+DPQ2+^-i 
= (47rr^Df)-Mn(r^/r<,) 



(6.3) 



leading to a logarithmic temperature dependence 
through T^p (T) . Similarly, contributions from the two di- 
agrams of type (b) can be evaluated to give 



-36 



^ jkajkl3fk.-k/af-k,klcr^ 



Hmp 



X {jkajk0fk,-k/crf-k,kla)k^ 
X {3ka3kl3fki,-kaf-ki,ka^^^ 

SO that 



(6.4) 



X {^kaJk'lAfk-klaf-k^klcr'^fkl-kaf-klM^''' ^^'^^ 

In a similar fashion, the total contributions from all dia- 
grams can then be written as 



a/3 



j: 



af3 



4j; 



a/3 



(6.6) 



where 



— ~2~0feaife/3) 

^Fa „ _ 

J2^ = {vW)-\fkJk'0[ft-kJ\,k,. 

+ fki,-kaf-ki,ka)k 

[/ k, — kiuf —ki,kfaf ki,ka 
~ f-ki,kufki,-kiafk,kia])k,ki,ki 

= {'fpl-y) \ijkajk'l3fk,k2<Tf-ki,kia 

X fk,,-k2afki,k(7)k,kf,ki,k2- {^■'^) 

Here J^^ corresponds to contribution from diagram (a) 
of Figure 3, J^^ is a sum of contributions from the two 
diagrams of type (b), J^^ is a sum of contributions from 
two diagrams of type (c) (the other two of type (c) gives 
J"^ = J^^) and J^^ is a contribution from diagram (d). 
In the above, we have used the relation {rUmp/T^N^) = 
l/(27ar<.). 



B. Phase relaxation rate 

The Cooperon contribution depends on the phase re- 
laxation rate r~^, which grows linearly with tempera- 
ture T. In general, this may be cut off by spin-flip scat- 
tering Tg, by spin-orbit scattering Tgo, or by a magnetic 
field characterized by ojh, ah of which are independent 
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of temperature. Therefore, a logarithmic temperature 
dependence in the conductivity requires that the phase 
relaxation rate satisfies the inequality 

max(l/rs, l/Tso,t^H) < I/t.^ < l/rtr- (6.8) 
The contribution to from e-e interaction is given by 



T epTtr 

In — - — 



<^FTtr 



(6.9) 



This is typically too small to satisfy the above inequality 
in thin ferromagnetic films where in particular the inter- 
nal magnetic field Bin can be estimated to give rise to 
ujh = '^{^FTtr)[^Bin / m* c) which can be large. A much 
larger contribution is obtained from scattering off spin- 
waves in such systems |15i] , which is given by 



1/r^ = 47rT 



J' 



(6.10) 



where J is the exchange energy of the s-electrons and 
Ag is the spin- wave gap. As estimated in Ref. [13], 
with this contribution to the phase relaxation rate, the 
inequality (6.8) can be satisfied within experimentally ac- 
cessible disorder and temperature ranges where the WL 
effects can be observed. 



VII. STRONG SHORT RANGE IMPURITY 
SCATTERING 

The results of the previous section can in principle be 
used to obtain the weak localization corrections to both 
longitudinal and Hall conductivities. However, The al- 
gebra gets fairly involved without contributing extra in- 
sight into the problem. Since higher angular momentum 
components are expected to be smaller, we will consider 
the dominant contribution that arises from a short range 
impurity model and show in the Appendix how effects of 
finite range anisotropic scattering can be included within 
model calculations. On the other hand, we will keep the 
calculations valid for arbitrary strength of the impurity 
scattering. 



A. Scattering amplitude, relaxation rate and 
particle-hole and particle-particle propagators 

These were already obtained for short range strong im- 
purity scatterings in Ref. [10.] and we will simply quote 
the results. The scattering amplitude is given by 



\/w Ju 



isujn + 2m(/c • kl)] . 

Here we defined w — u)/(l + u'), and u = uj (1 + w), where 
w = {■kN^VY and u = {g„/2Yw , and all quantities 
depend on the spin orientation a (suppressed here and 



in the following, except in the final expressions involving 
spin summation). In terms of the angular momentum 
components of / defined in (3.9), f^, we have from (7.1): 

/4, = /o + .f!k+kL + .fl,k-% /o - ^ - 

Jw 



/ii 



-isu ± 



= ; /;„ = 0, ]ml>l. 



(7.2) 



Using Eq.(7.2), the single particle relaxation rate given 
by Eq. (3.10) becomes 



2^ ~ 



(w + 2u). 



(7.3) 



One observes that is proportional to the Fermi en- 
ergy, the average number of impurities per electron and 
the dimensionless factor {w + 2u), expressing the effec- 
tive scattering strength per impurity. Eigenvalues of the 
particle-hole scattering amplitude t^jJ7 are obtained to be 

Aq = 1; A_frj = A„, 

Ai = 2wu{w + 2u)-\l + is—T^^) 

\2 = —{w + 2u)-^{u-1 + 2is^t';^) (7.4) 
u 

while for t^j^, one obtains (with = 

,1-tt. 



Co = {w + 2u)- 



'l + w 

^ -2wu{w + 2u)^'^{l + is—=) 



(1 — w — 2is\/w) + 2u 



w 



6 



-u{w + 2u) 



(7.5) 



It may be shown that A tkkiil) defined in (3.13) gives 
rise to small corrections to the diffusion coefficient, of 
order {I/eft) and hence may be dropped. 

Eigenvalues of the particle-particle scattering ampli- 
tude i^Jfc^" are obtained to be 



\l = [w-2u{l-2u)\/{w + 2u) 
\l, = [2wu±2^^r^^^)/{w + 2u) 



Aj.2 = ul{w + 2u) 



(7.6) 



We observe that Aq 1 if skew scattering is present, as 
it violates time reversal symmetry. 

The phase relaxation rate (t^°)^^ defined in Eq. (3.42) 
is given by 

(r^°)-i = r-i4u(l - u)/[w - 2m(1 - 2m)] (7.7) 

which is positive for not too large spin-orbit scattering, 
u < w/2, or (7cr < 1. 



B. Hall conductivity 

The conductivity tensor due to skew scattering was 
already evaluated in section IV. A for general strong finite 
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range impurity scattering in terms of the eigenvalues of 
the particle-hole propagator A. In particular, it gives 



A'/ 



'^r. 1 - A'l 



{7- 



For short range scattering, Eq. (7.4) gives explicit ex- 
pressions for the eigenvalues in terms of the scattering 
potentials. The side jump contribution was already eval- 
uated in Ref. [l3 and we quote the result: 



ZTT ^ — ' W -\- 2.U 1 + U 

a 

Using (7.4), this yields, in the small it ^ it; <C 1 limit 

1 



w (1 + A'l) 



(7.9) 



e2 ^ 



■^v 2ti ^ 1 - A' 



Weak localization correction 



(7.10) 



Evaluation of 3°"^ defined in section VI (Eqs. (6.6), 
(6.7)) in the present short-range (but arbitrary scattering 
strength) model gives 



-tXX 
Jy — 
TXX 

'H — 


(l + A'i)2-(A'i')2; Jf'' = 2A'i'(H- 
[A'l Jr - A'/Jf'']; - [A'/Jr + 


A'l) 

A'lJr^] 


-TXX 

J3 — 


^{22A'i Jr - (2m + 1)A'/Ji"^} 




i-xy _ 


^{(2m+ 1)A'/Jf" + 2SA'iJf^} 




TXX 

"^5 — 


-i{(2u- l)A'iJr - 2MA'/Jf^} 




i^y _ 

•^b — 


-i{2uA'i'jr + (2w- l)A'iJf^} 




^Jf - 


tq/3 1 7Q/3. 7Q/3 lap tq/3 


. (7.11 



VIII. COMPARISON WITH EXPERIMENTS 

Experiments measure the longitudinal and Hall resis- 
tances Rap as functions of both sheet resistance and tem- 
perature. In order to compare, we obtain the normalized 
relative conductances defined as 



1 Sa. 



Q/3 



Q/3 



(8.1) 



where Lqq = e^/27r^ and _Ro = ^/o'xx As shown above, a 
logarithmic temperature dependence in these quantities 
can arise either from interaction corrections or from weak 
localization corrections. However, although two separate 
groups have seen such logarithmic temperature depen- 
dences [1^ [l3| , the prefactors seem to be more universal 
for A^axx, independent of sheet resistance Rq or sample 
preparation for a range of Rq, but clearly disorder and 
sample dependent for A^cTxy in the same range of Rq. 
In this section we collect all our results above to obtain 
the total contribution to A^aap from all possible mecha- 
nisms considered above. As used in the text, superscripts 
ss and sj will refer to the skew scattering and side jump 
mechanisms, and / and WL will refer to the interaction 
and weak localization corrections, respectively. While the 
results for cr^^ and Sa^y are valid for finite range strong 
impurity scatterings, others are evaluated within a short 
range strong impurity scattering model. We have also 
assumed that the spin-orbit coupling is weak. 

The conductivities due to skew and side jump scatter- 
ings are 



A'l' 



^1-A'i 



e^^ (1^ 
2tt^ '^'^^''11- All 2 



■2) 



We may combine this into the compact expression 
J^^ = i?e{A}; ry^Im{A}; A ^ 



1-Ai 



(7.12) 



Note that the final result for J"^ contains detailed ef- 
fects of the potentials only through the eigenvalues A. 
This suggests that the results may be more general than 
the short range potentials used in the calculations. Also, 
as we will show in the Appendix, A'^ may approach unity 
in the limit of extreme forward scattering. 

In any case, for the short range impurity scattering 
model considered above, we then have contributions from 
weak localization corrections given by 



6a 



WL 

xy 



47r2 

Im{A) A'l' 



Re{A) 1 - A'l ' 



(7.13) 



Quantum corrections to the conductivities due to 
Coulomb interaction and weak localization effects leading 
to a logarithmic temperature dependence are 

fcr,'' = Loohxx HTt); Sal'x'^'^ = Loo HTr) 

r ssj n r ss.WL r ss.WL ^1 

da^y = U, da^y = da^^ ^_ y 

^ ""xy 



(8.3) 



The total conductivities and quantum corrections are 
simply 



ss . ss I SJ 

"xxi ^xy — a^y -r O^y 



Saxx = 5ar/ + <^a-^^^ <5a,, = Ja,^^ (8.4) 
Using these results, we obtain 



^xx 0(^xx + ^Oxx 
L\ CJxx — T 



(1 + hxx) In(rr) 
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A ax 



'xx " xy 

^00 a^'y + axi 



1 



(1 + r^y) 



In(Tr) (8.5) 



where h^^ defined in Eq, (5.23) is the exchange plus 
Hartree interaction contribution to the longitudinal con- 
ductivity and we have defined 



_ '^xy 



I xy 



(8.6) 



xy 



as the ratio of side jump to skew scattering contributions 
to the Hall conductivity. Note that Txy is a non-universal 
quantity. As shown in [14] , all current experiments can be 
understood if h^x <C 1 and r^y is sample dependent and 
is allowed to vary with disorder. In particular, this means 
that while the skew scattering and side jump mechanisms 
both contribute to the AH conductivity, the side jump 
contributions to the longitudinal conductivity as well as 
to the weak localization corrections to the conductivity 
tensor are much smaller than the corresponding skew 
scattering contributions when the spin-orbit coupling is 
weak. 



IX. SUMMARY AND CONCLUSION 

We develop a systematic general formulation for the 
AHE for strong, finite range impurity scattering start- 
ing from a microscopic model of electrons in a random 
potential of impurities including spin-orbit coupling. In 
particular, we consider quantum corrections to the AH 
conductivity, observed in different experiments on dis- 
ordered thin ferromagnetic films with apparently differ- 
ent results. General symmetry arguments presented here 
show that the e-e interaction corrections must vanish ex- 
actly, which then implies that there must be weak local- 
ization corrections in these ferromagnetic films despite 
the presence of large internal magnetic fields. 

Our evaluations of the WL effects within a short range 
but strong impurity scattering lead to the normalized rel- 
ative conductances given by Eq. (8.5), where the spin- 
orbit coupling has been assumed to be weak. These re- 
sults are consistent with all experimental observations, 
where the difference between different experiments arise 
due to different contributions from skew scattering vs side 
jump mechanism. 

In this paper we have only briefly mentioned the Berry 
phase effects. A systematic study of the Berry phase 
contributions to the AHE will be reported elsewhere. 
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APPENDIX: 



LONG RANGE CORRELATED POTENTIALS 

For completeness, here we consider models to incorpo- 
rate possible effects of small and large angle scattering. 



A. Model of small angle scattering 

Long range correlated potentials will scatter electrons 
predominantly by a small angle 9 « tt . A simple model 
is provided by a gaussian dependence 



v{\i - k/) = v{e) = 4V^yo(?(7'e-(«/«°)' 



(9.1) 



where Oq << n. The angular momentum components of 
V{9) are given by 



pTT in 

Jq 27r 



In the limit of weak scattering we have /, 
then 



o/4 (9.2) 

ma — 



{TTN„VofV2^/eo. 



(9.3) 



Neglecting skew scattering for the moment wc find 



m 



{m^ + {m-lY) 



(9.4) 



It follows that 1 — t^^ 
diffusion coefficient is enhanced by a factor 



o/8 <C 1 and therefore the 



DID, = {elis)-\ 



Model of strong back-scattering 



(9.5) 



It is well known, that the scattering of conduction elec- 
trons in amorphous metals can be anomalous in the sense 
that the transport relaxation time is smaller than the 
single particle relaxation time. This is due to the fact 
that the atomic structure is characterized by finite range 
order. The pair correlation function shows enhanced 
peaks corresponding to the nearest neighbor, next near- 
est neighbor, etc. shell. In other words, the system shows 
crystalline order over a certain usually short distance. As 
a consequence electrons are suffering Bragg scattering by 
large angles. The scattering cross section for large angles 
is larger than that for small angles. Consequently the 
angular average of the cross section <t{9), weighted with 
the factor (1 — cos 6*), appearing in the expression for 
the transport relaxation rate is larger than the uniform 
average in the single particle transport rate. In the case 
of polycrystalline material we expect a similar effect. 

The scattering potential V{r) of a crystallite or a small 
grain of amorphous metal will show oscillating behavior 
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in real space reflecting the nearly regular arrangement of 
atoms, and its Fourier transform will show a peak at a 
finite momentum q = lix j a corresponding to the spatial 
period a , which will be equal or close to the lattice con- 
stant of the crystalline phase. The width of the peak will 
be determined by the range of the short range order or 
the size of the crystallites. This is in contrast to a usual 
impurity potential whose Fourier transform has a peak at 
q — Q and a width corresponding to the range of the po- 
tential. In terms of the angular momentum components 
Vi of the scattering potential a peak in V^q) implies that 
some of the Vi will be negative. In particular, the com- 
ponent Aiof the t-matrix t^k, determining the transport 
relaxation rate will be negative. 

Let us consider a simple model of a crystallite of size 
L. Its scattering potential seen by a conduction electron 
of the matrix (assumed to be isotropic, as appropriate 
for an amorphous system) is something like 



Vi(x) = Focos(27rx/a)6'(L/2- |x| 
= VqSi{x)- Id 
y2(a;,t/) = V^S^{x)S^{y); 2d 

The Fourier transform of S\ {x) is given by 



(9.6) 



L \Kco&{K) sin(K) - Ksin(ii') cos(k)] 
'5i(^)-2 ^23^^2^ ' (9-7) 

where K = kLjl^ n = ixLja. Si{k) increases linearly 
with k at small k, has maximum at /e « 27r/a and de- 
creases as 1/k for large k. We may model this behavior 
by 



V2{k)=Vo 



kkn 



^2 _ ' 



(9.8) 



where fco = 27r/a. Using the relation of the transferred 
momentum k = k/ — kj to the scattering angle <f), k^ = 
2fc|-(l — cos(f)), where |k/_i| = kp, we get 



7] + COS ( 



(9.9) 



where V = Vo/{kFV2), V = kl/2kl - 1. 

The angular momentum components Vi may be calcu- 
lated as 



Vi 



In particular we find 
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COS(I0)F2(</'). 



(9.10) 



Vo = -y(?7- l)-^/^arctanJ^— > 

TT \ ?] — 1 

2— ji 

Vi = —V{ , arctan < 

TT V^- 1 

< (9.11) 



A + V2} 
r]-l 



In the limit t] ^ 1 the ratio of the / = 1 and 1 = 
components is given by Vi/Vq = —r]. We may estimate r] 

by assuming Z electrons in a unit cell of area resulting 
27rZ/a2 and therefore rj = 27r/Z-l. For Z « 2, 5 



in kp 



appropriate for a mixture of Fe'^'^ and Fe^^ one finds 
rj ~ 1.5 and then Vi/Vq « —0.6. In the following we 
will take the Vi to be given parameters, which may be 
negative. 

In order to keep the calculation simple we will neglect 
all angular momentum components with |Z| > 2. Defin- 
ing dimensionless quantities Vi = irN^Vi as before, the 
dimensionless scattering amplitudes are given by 

/J = T/o/(l + isVb); /ii,, = V±l,a/(l+?sV^±l,a) 
y±i,a = Vi±V^r^^. (9.12) 

Assuming weak spin-orbit scattering we may expand in 



f±l,a — Y 



+ isVx 



± (1 + isVtf^T^^ 



(9.13) 



The normalization factor 70 entering the expression for 
the relaxation rate is obtained as 

70 = T^ + 2-^+O(^), (9.14) 
l + w 1 -I- Wl 

- 2 - 2 

where w = Vq , wi = Vi . The eigenvalue Ai of tkki is 
found as 

^ ^ 1 21/0^1(1 + ^0^1) 
^ 70'- (l-hi«)(l -hwi) 

Analyzing this expression one finds that the largest nega- 
tive values of Ai are reached for weak scattering, Vo,Vi ^ 
1, when 



Ai 



2Vo 
w + Wl 



[yi + i(yo-2yi)\/^r^,] (9.16) 



The minimum of A'l is obtained if Vi / Vb = — 1/ \/2, where 
Ai =-l/V2. 

Let us now consider diagram w^, which is determined 



by the parameter given by 



JXX 

•^2 



7o-H(l+Ai)^(/o+/:i,.+/o-/ri,J+c.c.} (9.17) 



^1 — /o /+l,o- + /o /-l,cr 
(1 -I- W)(l -I- Wl) 

In the weak scattering limit, we have 



/3i = 61/70 



2Vn 



W + Wi 



[Vr-i{2Vi + Vo)V^T^J, (9.19) 



which differs from Ai only by the sign of the term Vq in 
the imaginary part, i.e. P[ = A^. 
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